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ABSTRACT 

We examine the relation between black hole accretion and bulge star formation as a 
fmiction of look-back time (r) in 20,541 obscured AGNs (with redshifts z ~0.1 and 
bolometric luminosities LboI ~ 10'''^-10'*^ erg s~^) optically selected from the Sloan 
Digital Sky Survey (SDSS). To quantify the most recently formed stars with ages less 
than typical AGN lifetimes, we estimate the differentiated specific star formation rate 
(SSFRt-) based on population synthesis analysis. Eddington ratio (A) is inferred using 
[O III] A5007 luminosity and stellar velocity dispersion as proxies for LboI and black 
hole mass respectively. We find that when t < tq, SDSS AGNs follow a power law 
A oc SSFR^°~^-^; the relation flattens out when r > tq. The threshold timescale tq is 
^ 0.1 (~ 1) Gyr in young (old) bulges. The scatter in the power laws is dominated by 
observational uncertainties. These results may provide useful constraints on models 
explaining the correlations between AGN activity and bulge star formation. 

Key words: galaxies: active — galaxies: evolution — galaxies: nuclei — galaxies: 
starburst — galaxies: stellar content. 



1 INTRODUCTION 

The growth of supemiassive black holes (SMBHs) could 
be linked with star formation (SF) in their host-galaxy 
bulges through multiple processes (see e.g.. iHec kman 2009|, 
for a review ) . Sta r s form under the cosm ological infall of 
gas (JBinnevI [l977i : iRees fc OstrikeJ [IqTtI '). some of which 
could be accreted by central SMBHs (e.g., iBurkert fc Silkl 
I2OOII ). Galaxy mergers or interactions may trigger starburst 
and quasar activity (e.g., (Sanders ct al. 1988; Moore et al. 
Il99a ). which could output energy and mom e ntum regulat- 
ing further growth (e.g., fsilk fc Ree^ Il998l ; iFabianI ll999l : 
ISpringel et al.l l2005f). AGN-driven outflows could boost 



starbursts (e.g., iBegelman fc Gioffil 1 19891 : ISilk fc NormanI 
I2OO9 ). After bei ng active, evolved stars may fuel AGNs via 
mass losses fe.g.. lNorman fc Scovillelll988l ; ICiotti fc Ostrikeil 
119971 ) . These processes may have different significance in reg- 
ulating SMBH growth and bulge SF in various AGN popu- 
lations at different cosmic epochs. 

While post-starburst populations have been ex- 

qua sars ( e.g., Boroson fc Okj 
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specific star formation rate (SSFR = SFR/Af,) and 
Eddington ratio (A = I/Boi/i/Edd). By studying the A 
distribution in AGNs as a fu nction of mean stellar age, 
iKauffmann fc HeckmanI l|2009l ) suggested that local AGNs 
occupy two distinct growth r egimes. These authors used 
the 4000- A break fD„(4000): [Bruzuall [1981 : iBalogh etTll 
ll999) to characterize mean stellar age. Arising from a 
series of metallic lines , Dn(4000) is most sensit ive to 1-2 
Gyr old stars (e.g., iBruzual fc Charloj 120031 ). whereas 
luminous AGN phase s appear to be brief (<0.1 Gyr; e.g., 
IYu fc Tremain3l2002h . 

To get higher temporal resolution and to quantify the 
most recently formed stars with ages less than typical AGN 
lifetimes, we examine SSFR averaged in the past r years 
(SSFElt) based on population synthesis analysis (e .g., iBical 
1 19881 : ICid Fernandes et alll2005l : lAsari et alllJOOTl ). Under 
the multiple-starburst approximation, the differentiated SF 
histories in a statistical sample of AGNs offered by the Sloan 
Digital Sky Survey (SDSS; York, ct al. 2000) enable us to in- 
vestigate the correlation between AGN and recent bulge SF 
activity as a function of look-back time. This approach offers 
insights otherwise unattainable with single and older age in- 
dicators. We describe the AGN sample and the assumptions 
in estimating A and SSFR^ along with uncertainties and se- 
lection biases in fj2] fj3] focuses on the correlations between A 
and SSFRt as a function of r. We discuss implications and 
caveats in ^ A cosmology with f7,„ = 0.3, Qa ~ 0.7, and 
/i = 0.7 is assumed throughout. 



X. LIU 



2 DATA AND METHODOLOGY 

2.1 The AGN Sample 

The AGNs in this study were drawn from the MPA 
SDSS-DR4 type 2 AG N catalog iKauffmann et al.1 1200 j 
iBrinchmann et alJlJOoi '). Type 2 AGNs are thought to be 
the obscured counterpa rts of type 1 AGNs according to 
unific ation models (e.g., I Antonuccil 1 19931 : lUrrv fc Padovanil 
Il99a ). The parent sample was selected from SDSS emission- 
line galaxies at z < 0.3 based on diagnostic line ratios that 
distinguish AGN from stellar ionizations (I Baldwin et al.l 
Ll981). The four relevant lines (H/3, [OlIl]A5007, Ha, 
and [N II] A6584) were detected with signal to noise ra- 
tios (S/N) > 3. We adopt the empirical criterion, 
log([Oiii] A5007/H/3)> 0.61/({log([Nii] A6584/Hq)-0.05)-F 
1.3), wh ich separates AGNs from star-forming galax- 
ies IjKau ffmann et al. 2003). The analysis is limited to 
[Qui] A5QQ7 luminosity I/[oiii] in the range of lO^'^-lO^'^Lo 
which corresponds to bolometric luminosities Lboi ^ lO'^^- 
10fl_erg s"£(e.g., iHeckman et~ai] l2004l : [Reyes et al.ll2008l : 
iLiu et al. 2009). The lower bound is to avoid faint-end in- 
completeness in a flux limited sample; It also effectively ex- 
cludes low ionization objects. The upper b ound is to avoid 
substantial scattered- AGN contamination (IZakamska et al.l 
l2006l : iLiu et al.l l2009l : iGreene et all |2009| ). which makes 
population analysis more complicated and uncertain. Re- 
sults on the L[oiii] > IO^'^Lq objects will be pre- 
sented elsewhere. The resulting sample consists of 20,541 
objects with a median redshift of 0.1. These include 
11,266 AGN-dominated objects with log([OlIl] A5007/H/3)> 
0.61/(|log([Nii] A6584/Ha) -0.47) + 1.19| abo ve the the- 
oretical "starburst limit" jKewlev et al.1 boOll ) and 9,275 
starburst-AGN composite^. 



2.2 Bolometric Luminosity and Black Hole Mass 

We use I/[oiii] (from the MPA catalog, uncorrected for ex- 
tinctiorjj) to infer intrinsic AGN luminosity. I/[oiii] is ob- 
served to be correlated with broad -band continuum lumi- 
nosities in unobs cured AGNs (e.g. iKaufFmann et al.l 120031 : 
iReves et al.ll2008l ). despite having a significant scatter (e.g., 
^0.36 dex in logL[oiii] in the correlation with M2500 
iReves et al.ll20 08V We adopt the conversion log(LBoi/i0) = 
0.99 X log(L[o HI]/!/©) + 3.5( ±0.5) appropria te for I/[oiii] 
before extinction correction (jLiu et al.l |2009| ). We correct 
L[oiii] for the contribution from star formation using an em- 
pirical method based on an object's positi on on the [N II]/Hq 
vs__[0lIl]/H/3 diagram following Kauffmann fc HeckmanI 
(|2009l ). SMBH mass is estimated (Mbh) from stellar velocity 
dispersion (<t*), a ssuming the relation observed in local inac- 
tive g alaxies fe.g.. lGebhardt et al.ll200d:lFerrarese fc Merritj 
l20od ) calibrated bv iTremaine et al.1 (|2002l '). We adopt ai 



^ It is important to include starburst-AGN composites in the 
analysis because (i) they are a significant population; and (ii) 
otherwise the sample would be biased against galaxies with high 
star formation and low AGN activity. 

■^ We adopt iroiiil uncorrected for extinction and the appro- 
priate calibration rather than the extinction-corrected method 
because the latter approach is sensitive to uncertainties in the 
assumed red dening laws and re lies on over-simplified dust-screen 
models (e.g.. lReves et al.ll200^ . 



from the SDSS specBS pipeline (|Adelman-McCarthv. et al.l 
[2OO8) and exclude uncertain a* measurements (< 10 or 
> 600 km s^M. 



2.3 Differentiated Specific Star Formation Rate 

The common emission-line SFR indicators cannot be di- 
rectly applied to luminous AGNs due to AGN contamina- 
tion. Continuum indices are adopted such as D„(4000) (e.g., 
[Kauffmann et al. 2003 : Brinchmann et al. 2004 ). To achieve 
higher temporal resolution than D„(4000), we estimate 
SSFR averaged in the past t years (SSFRr) from popula- 
tion s ynthesis analysis (e.g.. lBicalll988l : ICid Fernandes et al.l 
[2005). Compared to SFR, SSFR is less subject to absolute 
calibration un certainties associated wit h fiber effects and ex- 
tinction (e.g.. IBrinchmann et al.ll2004l ). We measure SSFRs 
inside the SDSS 3"fibers (corresponding to the central ~1.5 
kpc at typical redshift of the sample) to quantify the bulge 
SF. 

We adopt continuum m odels from the MPA DR4 
archive IjTremonti et al.ll2004l ). and have t ested them on a 
randomly drawn subsample with our code (ILiu et al.ll2009l ). 
The assumption is that the bulk stellar population was 
built up in multiple starburst events with ages smaller than 
the age of the universe at the galaxy redshift. To sam- 
ple the whole history of a host galaxy, stellar continuum 
is fit as a linear combination of ten instantaneous star- 
burst models (broadened with the measured ct,) with ages 
of 0. 005, 0.025, 0.10, 0.29, 0.6 4, 0.90, 1.4, 2.5, 5.0, and 11 
Gyr (jBruzual fc Charlotll2003l '). along with dust attenuation 
as an additional parameter. There are three sets of mod- 
els with different metallicities (0.2, 1.0, and 2.5 solar). We 
adopt the metallicity with the smallest x^ as the best es- 
timate. Given a best-fit model, SSFR-r is calculated as the 
summed ratio of the mass and age of the starburst templates 
that fall in interval r, divided by the total mass contained 
in all templates. The model coefficients in the MPA DR4 
archive are in terms of luminositj|j, and need to be con- 
verted to mass weight s using the mass-to- l ight r atios of the 
starburst models from lBruzual fc CharlotI l|2003h . 

We study the relation between AGN activity and 
recent bulge SF over three characteristic time inter- 
vals. First we divide the sample into young and old 
bulges, because the distribution of m ean stellar age is bi- 
modal l|Kauffmann fc HeckmanI l2009l i. For "young" bulges 
(D„(4000) < 1.6JJ, the characteristic times are r < 5 
Myr (containing the age grid 0.005 Gyr), ~ 5-100 Myr 
(0.025 and 0.10 Gyr), and ~ O.I-l Gyr (0.29, 0.64, and 
0.90 Gyr). The large separations (~l-dex) between these 
time intervals ensure the recovered SF at different e pochs to 
be w ell distinguished from one another (e.g., lOcvirk et al.l 
I2OO61 ). Higher temporal resolution would cause over inter- 
pretation due to degeneracies between models with sim- 
ilar ages. The timescales are characteristic of: 1-10 Myr 



•^ We caution that the SSFRs in IChen et al.l l l 20091) appear to 
be inferred using the luminosity coefficients directly without the 
requisite mass-to-light conversion. 

* Our results are independent of the specific D„(4000) division 
value chosen in the range of 1.6 to 1.8, but we adopt 1.6 so that 
the old bulges contain sufficient objects for reliable statistics. 
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Figure 1. Diflferentiated specific star formation rate (SSFR^) versus Eddington ratio A. Tlie panels (a)-(c) are for young bulges (Dn(4000) 
< 1.6) whereas tiie panels (d)-(f) are for old bulges (D„(4000) > 1.6). Look-back time intervals r are labelled on each plot. Data and 
distribution curves are color coded in SSFR,- (also discriminated with line thickness). The numbers (also color-coded in SSFR,-) are the 
Kolmogorov-Smirnov (KS) significance denoting the probability that a distribution is normal. Red contours represent ten evenly spaced 
levels in number density. Black solid and dotted lines show the bisector linear regression fits (median and 1-cr uncertainties; Table [T}. 



for massive star lifetimes (with ma sses > 10 Mq), <0.1 
Gyr for luminous AGN lifet ime (e.g.. lYu fc Tremaind[20o3 : 
iHopkins fc Hernguistl 120091 ), and 1-2 Gyr for lifetime of 
lower-luminosity AGNs (e.g.. lHopkins fc Hernauistll2009l) to 
which Dn(4000) is m ost sensitive to (e.g.. lBruzual fc CharlotI 
l2003l : iKauffmann e t al. 200i). For "old" bulges (D„(4000) 
> 1.6), there is little current SF, and SSFRr cannot be re- 
liably estimated for r <g 0.1 Gyr (little mass allowed in the 
youngest starburst grids). The adopted time intervals are 
r < 0.1 Gyr (containing the age grids 0.005, 0.025, and 0.10 
Gyr), ~ 0.1-1 Gyr (0.29, 0.64, and 0.90 Gyr), and 1-2 Gyr 
(1.4 Gyr). 

2.4 Uncertainties and Selection Biases 

The uncertainty of A is dominated by the systematics of 
LboI and A/bh. The quoted 0.5-dex I-ct uncertainty on LboI 



contains the intrinsic scatter of the L 



[OIII]~-^2500 



relation 



and th e uncertainty of bolometric correction in unobscured 
AGNs (jLiu et al.ll2009l ). For Mbh, the uncertainty is propa- 
gated from at added in quadrature with the 0.3-dex scatter 
of the Mbh-o- relation l|Tremaine et al.ll2002l ). The total 1- 
a uncertainty of A is ~ 0.6 dex. We have tested effects of 



Malmquist biases due to these scatters, and found that there 
is little bias on Mbh, since the a^ distribution is close to log 
normal. On the other hand, there is a small yet detectable 
bias on LboI, such that the inferred LboI is an overestimate 
(by < 0.03 dex) of the true value (assuming LboI is a prior 
whereas Z/[oiii] is a posterior) gi ven an observed bottom- 
heavy distribution of I/[oni] (e.g.. lShen et al.ll2008l ). 

We have empirically tested the continuum-SSFRr ap- 
proach using 12,090 S PSS star- forming galaxies (e.g., 
iBrinchmann et al.ll2004l ). Our initial calibrations were car- 
ried out on high-S/N objects (with median S/N > 20 
pixel" ^ and S/N > 5 in a*). To constrain selection bi- 
ases we tested low-S/N objects using co-added spectra 
and found similar results. We have verified that SSFRr 
correlate with the current SSFRe usin g the MPA SFRe 
estim ates from emission-line modeling (|Brinchmann et al.l 
l2004h . SSFR<5Myr ~ SSFRe with a 1-cr scatter of -0.25 
dex in a bisector fit; SSFRr with larger r correlates with 
but does not equal to the current SSFRe, probably indicat- 
ing varying SF as a function of look-back time. While these 
tests are limited to star-forming galaxies, they do lend sup- 
port to our general approach, especially for young bulges. 
We take the 0.25-dex scatter as our fiducial estimate for 
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the SSFRr systematic unc ertainty. Th i s is c omparable to 
the 0.3-dex estimate from lAsari et al.1 (|2007l l by compar- 
ing continuum-based SFRs against Ha SFRs for SDSS star- 
forming galaxies. We have added noise to spectra of a ran- 
domly selected subsample, and found that typical errors of 
SSFRs due to noise in the continuum spectra are no larger 
than ~ 0.2 dex. We take 0.3 dex as our fiducial estimate for 
the total l-a uncertainty for SSFRr . 



3 RESULTS 

Fig. [T] shows Eddington ratio (A) versus SSFR-r, and the A 
distribution (P(logA)) as a function of SSFRt . The young 
bulges (D„(4000) < 1.6) have median log(I/[oiii] / Lq) = 
6.89 (with a I-ct scatter of 0.39 dex) and log(M./M0) = 
10.88 (±0.32), whereas the old bulges (D„(4000) > 1.6) have 
median log(L[oiii] / Lq) = 6.73 (±0.35) and log(M*/M0) = 
11.13 (±0.33). For t ^ to {tq ~ 0.1 Gyr for young and ~ 1 
Gyr for old bulges), both populations appear to follow power 
laws. The correlations are strong: the Spearman probabili- 
ties for null correlation are close to zero (Table 1, Col. 8); 
the relations flatten out when r > ro (Table 1, Col. 7). The 
cutoff in Fig. [if at the high SSFRr end is artificial due to its 
definition and the limited temporal resolution of the contin- 
uum approaclij. 



3.1 The A-SSFR Relation 

A of AGNs in young bulges correlate with their current 
SSFRs (SSFRr<5Myr; Fig. [TJi). We fit the relation A = 
10''(SSFR/Gyr^^)'°, the linear regression results of which 
are given in Table \T\ in. We adopt bisector fits appropriate 
for two variables with no assumption of one depending on 
another. The fit for young bulges is logA = (—0.42 ± 0.01) ± 
(1.12±0.01)log(SSFR/Gyr~^). The scatter (1-cr orthogonal 
of 0.54 dex) is dominated by observational uncertainties. 
Given any fixed SSFRt, P(A) is roughly log normal, and its 
central value increases with increasing SSFR^. While some 
of the Kolmogorov-Smirnov (KS) probabilities labelled in 
FigH] are small, when we divide the sample into finer bins 
in SSFRt the KS probabilities become larger. 

In old bulges with measurable SSFRs, A also correlates 
with recent SSFRs (SSFRT<o.iGyr; Fig.[T]i). The fit for old 
bulges is logA = (0.00±0.01) + (1.00±0.01)log(SSFR/Gyr-^) 
(Table [l]). Considering uncertainties, this is similar to the 
relation in young bulges, although the median A is ~ 1 
dex smaller than in young bulges. Similarly, given any fixed 
SSFRt, P(A) is close to log normal, and the central A in- 
creases with increasing SSFRt. 



^ The cutoff is reached if all mass is recovered in the 1.4 Gyr 
grid. In practice this only happens for t >1 Gyr. For younger 
ages, the cutoff is never reached, because while most of the light 
may go into young populations, it is always the old populations 
which dominates the total mass (due to the drastically different 
mass-to- light ratios of old and young stars). 



Table 1. Results of the fits log(A) = b + k log(SSFRT/Gyr-l) 



Sample 


r/Gyr 


b 


k 


(7 


N 


P 


-Pnull 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


young 


^ 0.005 


-0.42 


1.12 


0.54 


9903 


0.23 


< 10-40 


young 


0.005-0.1 


-0.45 


1.03 


0.63 


8420 


0.15 


< 10-40 


young 


0.1-1 


-0.81 


1.00 


0.68 


9951 


0.01 


0.5 


old 


^0.1 


0.00 


1.00 


0.52 


3921 


0.29 


< 10-40 


old 


0.1-1 


-0.82 


0.78 


0.54 


4302 


0.29 


< 10-40 


old 


l-2t 


-0.37 


1.04 


0.48 


5351 


0.16 


10-16 



Col. (1): D„(4000) < 1.6 for "young" whereas D„(4000) > 1.6 for "old" 
bulges; Col. (2): look-back time t; Cols. (3) & (4): normalization and slope 
from the bisector fits with 1-a uncertainties of <0.01 dex; Col. (5): standard 
deviation of the distribution, a is measured along the orthogonal direction 
to the fit; Col. (6): number of objects in each subsample. To be included, a 
galaxy needs to contain mass in at least one of the starburst grids in interval 
t; Cols. (7) & (8): Spearman correlation test coefficients. | for Fig.[TJ, the fits 
are for 2499 objects with log(SSFR/Gyr— i) < —1, to minimize effects from 
the artificial cutoff at the high SSFR end. 

3.2 Correlation between AGN Activity and 
Recent Bulge Star Formation 

In young bulges A also correlates with SSFRs averaged over 
the past ~ 5-100 Myr (Fig. [TJj). The correlation for r ~ 5- 
100 Myr is similar to r ~ 5 Myr (Table[l]). For young bulges, 
the correlation flattens out when r ~ 0.1-1 Gyr (Fig. [l};): 
given a flxed SSFRt, P(A) is still largely log normal, yet 
A depends little on SSFRt averaged over the past ~ 0.1-1 
Gyr. Figs, [if e,f ) show similar results in old bulges with the 
main difference of having a larger tq (~ 1 Gyr). While there 
is still some correlation for r ~ 1-2 Gyr (Figs, [l}') partially 
due to the artificial cutoff, the relation starts to fiatten out 
(Spearman p being smaller; Table [TJ. 



4 DISCUSSION 

Our results imply similar relations between AGN activity 
and recent SF in old and young bulges, despite different Ed- 
dington ratios observed. The current Eddington ratio corre- 
lates with the most recent bulge SF: A oc SSFRt estimated 
in a certain look-back time interval r; it depends little on 
SF happened before some threshold time ro. The main dif- 
ference between young and old bulges is on ro : young bulges 
have ~ 0.1 G yr whereas old bulges have ~ 1 Gyr. 

Recently iKauffmann fc HeckmanI (|2009l ) have studied 
the A distribution as a function of D„(4000). These authors 
find that P(A) is log normal in young bulges and exponential 
in old bulges, and sug gest two growth regimes. W hile this 
could still be the case, JHopkins fc Hernguistl l|2009t ) suggest 
that different A distributions may also arise from universal 
AGN light curves with different quench times: AGNs in old 
bulges have apparent exponential P(A) because accretion 
and star forma tion were qu enched in a series of earlier epochs 
(see also e.g., IShenll2009l l. Our results point to the latter 
explanation. Namely, the physical processes that link AGN 
activity and bulge SF seem to work universally in young and 
old bulges, albeit with different quenching/threshold times. 

Using population synthesis analysis, we have found that 
when grouped in terms of similar SF activity with ages less 
than typical AGN lifetimes, BHs in old bulges also have log 



BULGE STAR FORMATION HISTORY IN AGNs 



normal P(A). We have shown that for BHs in young bulges, A 
does cor relate with the most recent SF (Figs.[l^,b). The rea- 
son whv lKauffmann fc HeckmanI l|200 9 I found little correla- 
tion is most likely that the temporal resolution of D,i(4000) 
is not enough to disintegrate the most recent SF in young 
bulges. The reciprocal SSFR^^, which is the mass- weighted 
stellar age averaged among the stellar populations formed 
over r, can be estimated for different timescales r, whereas 
D„(4000) is most sensitive to r = 1-2 Gyr. If we assume 
the observed sample consists of AGNs at different evolution- 
ary stages, the A-SSFR results (Table [l| imply light curves 
A oc i~' ■ ~ '^ , which is noticeably less steep than predic- 
tions A oc ^-(i-5-2 0) f]-o)2i niodels assuming ex treme gas ex- 
pel by AGNs fe.g.. lHopkins fc Herna"uistll2009l '). This milder 
decay could suggest more moderate self-regulated AGN and 
SF activity, although the observational uncertainties in our 
results are significant. 

Another implication is that AGN environmental stud- 
ies should draw control samples matched not only in the 
mean stellar age (as usually implemented using D„(4000)), 
but also in the most recent bulge SF properties. The effect 
would be particularly important for AGNs in young bulges, 
as we have shown that their SF-AGN-correlation timescale 
is usually much shorter (~ 0.1 Gyr) than that probed by the 
overall mean stellar age. 

Our results are limited to optical AGNs with z ~0.1 and 
i/Boi ~ lO^'^-lO*^ erg s~^, and may not apply to other AGN 
populations. In principle the approach could be extended to 
lower-luminosity AGNs (LINERs) and to higher-luminosity 
AGNs and quasars, although selection incompleteness and 
biases may be more severe. Multi-wavelength SFR indica- 
tors that probe different SF timescales or phases may also 
help elucidate the AGN-SF relation as a function of look- 
back time. The approach demonstrated based on population 
synthesis analysis relies on statistical AGN samples with 
high-S/N spectra. We look forward to the next generation 
of ground-based near-IR multi-object spectrographs (e.g., 
iBell. et al.ll20ld ) that will make an equivalent study possible 
at high redshifts to better understand the coupled growths 
of stellar bulges and SMBHs in the early universe. 
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